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The nature of superconducting fluctuation effects in the isovalently-substituted iron pnictide 
BaFe 2 (Asr-a;Pa ;)2 {x « 0.35) is probed through measurements of the magnetization and magneto¬ 
conductivity around the superconducting transition. The results, obtained with magnetic fields up 
to 9 T applied in the two main crystal directions, are consistent with anisotropic Ginzburg-Landau 
(GL) approaches for finite applied magnetic fields. The analysis allowed to determine with accuracy 
the out-of-plane, ^c(O), and in-plane, Ca6(0), GL coherence lengths. Significant differences are found 
between the ^c(O) values resulting from electrical transport and magnetization data. According to 
recent theoretical approaches, these differences could be interpreted in terms of the multiband na¬ 
ture of this material. The analysis of data in the low held region around the transition temperature 
also suggests that phase fluctuations, although possibly relevant in other Fe-based superconductors, 
may play a negligible role in this compound. 

I. INTRODUCTION 

Superconductivity appears in iron pnictides when the 
antiferromagnetic phase is suppressed by chemical dop¬ 
ing of the parent compound or by applying pressure 
It is now well known that an isovalent chemical dilu¬ 
tion (e.g., the partial replacement of As by P, or Fe by 
Ru)^A presents notable differences with charge doping 
(e.g., by partially replacing O by F in ROFeAs, where R 
is a rare earth element^^— or Ba by K in BaFe 2 As 2 )^. 

While isovalently-diluted compounds are closer to the 
clean limit charged dopants create strong scattering 
potentials that affect superconducting properties like 
the vortex pinning^i^ the upper critical field^^A^ and 
even the superconducting gap symmetry] ^^ d 
many studies indicate that iron-pnictides with charged 
dopants present a fully gapped Fermi surface (see e.g., 

SmFeAsOi-xFa;^ PrFeAsOi-j;^, BaFe2-xCoa;As2)^ 

Bai_xKa;Fe 2 As 2 )^ and BaFe 2 -xNia;As 2 ^i^), while iso- 
valently substituted compounds like Ba(Fei_xRua;) 2 As 2 
and BaFe 2 (Asi_a;Pa :)2 present evidences of nodesi^^— It 
is expected that fluctuation effects around the supercon¬ 
ducting transition may also present differences in these 
two types of iron pnictidesi^ These effects have been al¬ 
ready investigated in presence of charge-doping^^— but 
to our knowledge still remain unexplored in isovalently 
substituted iron-pnictides. 

The main aim of the present work is twofold: on the 
one hand, we will probe the nature of superconducting 
fluctuation effects in isovalently-diluted compounds, and 
the applicability of the existing Ginzburg-Landau (GL) 
approaches in the different regions of the H — T phase 
diagram. On the other hand, these effects will be used 
to obtain precise information about superconducting pa¬ 
rameters like the coherence lengths (or the upper critical 


fields), the anisotropy factors, and the system dimension¬ 
ality (a technique sometimes referred to as fluctuation 
spectroscopy^) . Fluctuation effects may also provide in¬ 
formation about the multiband nature of these materials. 
In fact, it has been shown that a two-band model pre¬ 
dicts a change in the relative amplitude of different fluc¬ 
tuation observables in relation with the single-band case, 
mainly through a renormalization of the c-axis coherence 
length— Among the isovalently-substituted iron pnic¬ 
tides, BaFe 2 (Asi_a;Pa ;)2 with x ss 0.35 presents the high¬ 
est critical temperature (Tc ^ 30 K) and thus is the best 
candidate to investigate superconducting fluctuations in 
these materials. In this work we present a systematic 
study of fluctuation effects around Tc in the electrical 
conductivity and magnetization of several high-quality 
single crystals of this compound. The measurements are 
performed in magnetic fields up to 9 T, which close to 
Tc is well inside the finite-field (or Prange) fluctuation 
regime— The results will allow to probe the applicabil¬ 
ity of GL approaches for finite applied magnetic fields, 
and to investigate possible differences in the coherence 
length amplitudes associated with the multiband nature 
of this material. Finally, the low-field (~ 10“^ T) behav¬ 
ior of the magnetization for temperatures around Tc will 
be studied to check the relevance of phase fluctuations in 
BaFe 2 (Asi_ 2 ;Pa;) 2 , which were claimed to be important 
in other iron pnictides— 


II. GROWTH AND CHARACTERIZATION OF 
THE CRYSTALS 

The BaFe 2 (Asi_xPx )2 single crystals used in this work 
were grown by using the Ba 2 As 3 / Ba 2 P 3 self-flux method 
described in Ref. Some details of their character- 
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FIG. 1. a) Example of a typical EDX spectrum, b) Example 
of x-ray diffraction pattern on a single crystal, obtained by 
using the geometry to observe the reflections by the ah layers. 
Inset: rocking curve associated to the (004) reflection, show¬ 
ing that the dispersion in the orientation of the crystal c-axis 
is about 0.1°. 


ization may be seen in Ref. Is^ . They are plate-like, 
with typical surfaces of several mm^ and thicknesses 
up to ^ 0.1 mm. Their stoichiometry was checked 
with a scanning electron microscope (Zeiss FESEM Ul¬ 
tra Plus) equipped with a EDX (Energy Dispersive x- 
ray) spectroscope. A typical EDX spectrum is shown 
in Fig. [iKa). The average stoichiometry resulted to be 
Bai.o 4 Fei. 9 iAsi, 33 Po. 72 , with a variation smaller than 
0.4% between the different crystals and the different stud¬ 
ied areas. The partial substitution of As by P is about 
35 %, which is close to the value that maximizes T^- 
The crystallographic structure was studied in some of the 
crystals by x-ray diffraction (XRD) by using a Rigaku 
MiniFlex II diffractometer with a Cu-target. A typical 
example of the reflections by the ab planes is presented in 
Fig- mb). The absence of reflections other than the (001) 


TABLE I. Physical parameters of the crystals. R1 and R2 
were used in the resistivity measurements, and the stack in 
the magnetization measurements. 


crystal 

surface 

(mm2) 

thickness 

(pm) 

mass 

(mg) 

R1 

0.78 X 0.49 

15 

0.036 

R2 

1.60 X 1.30 

18 

0.234 

Stack 

~5.20 

879 

29.514 


indicates an excellent structural quality of the crystals. 
An example of the rocking curve for the (004) reflection 
is presented in the inset of Fig. mb)- It confirms that 
the crystal c axis presents a dispersion of only ^ 0.11°. 
Powder x-ray diffraction in some grounded crystals al¬ 
lowed to determine the lattice constants of the tetragonal 
structure, which resulted to be a = 5 = 0.39249(15) nm 
and c = 1.2833(4) nm, in agreement with data in the 
literaturej^i^i^ 


III. RESISTIVITY MEASUREMENTS 

The in-plane resistivity (along the ab layers), poh, was 
measured in two single crystals with a Quantum Design’s 
Physical Property Measurement System (PPMS), by us¬ 
ing four contacts with an in-line configuration and an ex¬ 
citation current of ^ 1 mA at 23 Hz. The measurements 
were performed with magnetic helds up to fJ-oH = 9 T 
applied both parallel and perpendicular to the ab layers. 
The size of the crystals chosen for these measurements 
are presented in Table HI The finite size of the electrical 
contacts leads to an uncertainty in the pab amplitude of 
~ 25%. 

An example (corresponding to crystal 2) of the 
Pab(T)H behavior around the superconducting transition 
is presented in Fig. [2l In the overview shown in the in¬ 
set it may be seen a linear temperature dependence up 
to 100 K, without the kinks associated to structural and 
magnetic transitions typical of underdoped samples;^ In 
absence of an applied field, the resistive transition in both 
crystals is very sharp: the transition midpoint is about 
~ 28 K with a full width of 0.8 K as estimated from 
the 90%-10% criterion. This would allow to investigate 
fluctuation effects in a wide temperature range above the 
transition. However, the resistivity rounding just above 
Tc typical of fluctuations is almost inappreciable (see the 
inset in Fig. mb)). This is consistent with the relatively 
large normal-state (or background) in-plane conductivity 
aj^bJ S'® compared with the fluctuation-induced conductiv¬ 
ity predicted by the Aslainazov-Larkin approach for 3D 
materials, 

fl ^ I'D 

3271 ^ 0 ( 0 ) ’ ^ ^ 

where e = ln(T /Tc) is the reduced temperature, e is 
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FIG. 2. Example (corresponding to crystal 2) of the tem¬ 
perature dependence of the in-plane resistivity around Tc, for 
magnetic fields up to 9 T applied in the two main crystal di¬ 
rections. Inset in (b): detail of the normal-state behavior for 
iil = 0 up to 50 K. The line is the background contribution 
obtained by a linear fit above 35 K. No appreciable rounding 
associated with fluctuations is observed above the resistive 
transition. 


upper critical field of the band with the largest coher¬ 
ence length is about ^ 1 T. Note, however, that in the 
general case the different bands cannot be directly as¬ 
sociated with different coherence lengthsj ^^i^° The solid 
lines in Fig. are linear fits for iJ > I T. The resulting 
Hc 2 {T) slopes and the extrapolated Tc values are sum¬ 
marized in Table HU From these values, the in-plane and 
out-of-plane GL coherence length amplitudes, ^ab(O) and 
^c(O) respectively, were obtained through 


and 


U( 0 ) = 


00 


2TTT,fXo\dH^2/dT\ 


( 2 ) 


?c( 0 ) = 

where 

dHj^/dT 
^ ~ dH^JdT 

is the superconducting anisotropy factor. The resulting 
values are also compiled in Table HIl They are expected 
to be inappreciably affected by the small distribution: 
i) The percolative nature of the resistivity measurements 
would increase the measured Tc value to be introduced 
in Eq. But this effect is expected to be 0.4 K (the 
transition half-width, see above), less than 2%. ii) The 
magnetic field displaces the p(T)h curves to lower tem¬ 
peratures without introducing an appreciably broadening 
(apart from the rounding observed in the upper part of 
the curves for H || c). Then, the Hc 2 {T) slopes deter¬ 
mined from the transition midpoints (to be introduced 
in Eqs. m and O) are expected to be unaffected by 
a Tc distribution. The linear Hc 2 iT) behavior obtained 
above I T for both field orientations confirms this point. 


( 3 ) 

( 4 ) 


the electron charge, U the reduced Planck constant, and 
^c(O) the c-axis coherence length amplitude. By using 
the ^c(O) = 1 nm (see below), at e = 0.1 one obtains 
a^b ~ 2 X 10^ (Hm)“^, which is 10^ times smaller than 
aj^b ~ 1-3 ^ 10® (f2m)“^ just above Tc. It is worth noting 
that under the largest H amplitudes a slight rounding 
may still be observed in the upper part of the Pab(T) 
curves, that may be attributed to critical fluctuations 
near the Hc 2 {T) line. 

The absence of important fluctuation effects in the re¬ 
sistivity allows to determine with accuracy the tempera¬ 
ture dependence of the upper critical fields from the tran¬ 
sition midpoints. The result for the two studied samples 
is presented in Fig. |3l For both field orientations Hc 2 {T) 
is linear to a very good approximation, although for fields 
below ^ 1 T a slight positive curvature is observed (see 
the inset in Fig. [SKb)). This effect occurs at tempera¬ 
tures within the transition width and could be due to a 
Tc distribution. According to Ref. [H it could also be at¬ 
tributed to the multiband nature of this material, if the 


IV. MAGNETIZATION MEASUREMENTS 

The magnetization measurements were performed with 
a Quantum Design’s SQUID magnetometer (model 
MPMS-XL) in a stack of ten high-quality single crys¬ 
tals glued with a minute amount of GE varnish, with a 
total mass of 29.514 mg and a volume of 4.71 mm® as 
determined from their theoretical density. A picture of 
the individual crystals used in the experiments and of 
the resulting pile may be seen in Fig. 2] As we will see 
below, such a large sample is necessary to attain the reso¬ 
lution needed to study the weak diamagnetism due to su¬ 
perconducting fluctuations. By using optical microscopy 
we estimated that the misalignment between the differ¬ 
ent crystals in the stack is less than 2 °. As the crystals 
anisotropy is moderate (7 w 1 . 8 , see Table HIl) the effect 
of such misalignment is negligible: the relevant parame¬ 
ter is the upper critical field, which according to the GL 
theory depends on the orientation (relative to the c axis) 
as Hc 2 {d) = T[^/\/ cos^ 9 7 “^ sin^ 0 . Then, the error 
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TABLE 11. Superconducting parameters of the crystals studied, indicating the observable used to obtain them. Note the 
difference between the ^c(O) values determined from the Lf-dependence of the resistive transition, and from the analysis of 
flnctuation effects in the magnetization. See the main text for details. 


Crystal 

Observable 

p.odHc2/dT 

godUljdT 

Tc 

Ca6(0) 

Cc(0) 

7 



(T/K) 

(T/K) 

(K) 

(nm) 

(nm) 


R1 

p(T)h 

-1.72(2) 

-3.84(6) 

28.5(1) 

2.59(2) 

1.16(4) 

2.23(6) 

R2 

p(T)h 

-1.72(2) 

-4.05(9) 

27.3(2) 

2.65(3) 

1.13(5) 

2.35(8) 

Stack 

Mfi{T)H (Gaussian region) 

-1.81(2) 

-3.16(12) 

28.1 

2.54(3) 

1.45(6) 

1.75(5) 

Stack 

Mfi{T)H (critical region) 

-1.75(15) 

-3.2(3) 

28.2(1) 

2.6(1) 

1.4(3) 

1.8(3) 



T (K) T (K) 


FIG. 3. Temperature dependence of the upper critical field for 
H A- ab and H || ab, and for the two crystals studied. These 
data were obtained from the midpoint of the corresponding 
resistive transitions. Inset in (b): detail of the non linear 
behavior below 1 T. 


in the iJc 2 value is as small as 0.04% (for H || c) and 
0.14% (for H II ab). Even in the case that the misalign¬ 
ment were as large as 10°, the errors would be only 1% 
(for H \\ c) and 3% (for H || ab). 

The temperature dependence of the zero-field-cooled 
(ZFC) magnetic susceptibility, measured with a low field 
(0.2 mT) perpendicular to the ab layers, is presented in 
Fig. [5] These data are corrected for demagnetizing effects 
by using the demagnetizing factor needed to attain the 
ideal value of-1 at low temperatures {D = 0.71), which is 
consistent with the physical dimensions of the pile. From 
these curves, Tc was estimated by a linear extrapolation 
to X = 0 of the higher-slope x(T) data, and the transition 
half-width as ATc = Tco — Tc, where Tco is the highest 
temperature at which a diamagnetic signal is resolved 
(the procedure is detailed in the inset in Fig.[5|). In spite 
of the large volume of the sample, the ATc value is only 
0.4 K (the full width being 0.8 K), confirming its excellent 
stoichiometry, and allowing to study fluctuation effects in 
a wide temperature region above Tc{H). 

To measure the weak magnetic moment due to super- 



FIG. 4. a) Single crystals used in the magnetization measure¬ 
ments. b) and c) Top view and, respectively, side view of the 
pile formed with the single crystals shown in (a). 


conducting fluctuations above Tc {m ^ —10“® emu, see 
below) we used the Reciprocating Sample Option (RSO). 
We averaged eight measurements consisting of 10 cy¬ 
cles at 1 Hz frequency, which lead to a resolution in 
the ~ 10“® emu range. The magnetic fields used in 
the experiments range from 1 to 6 T which allowed to 
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FIG. 5. Temperature dependence of the low field (0.2 mT) 
magnetic susceptibility of the pile of single crystals. These 
measurements were performed with H ab after zero-field 
cooling (ZFC, circles) and field-cooling (FC, triangles). The 
data are already corrected for demagnetizing effects. Inset: 
detail around Tc showing the sharp diamagnetic transition 
(only 0.8 K full width). See the main text for details. 


deeply penetrate in the finite-field (or Prange) fluctua¬ 
tion regime (see below). The m[T) data around for all 
fields amplitudes and orientations studied are presented 
in Figs.|6l(a) and (b) (already subtracted from the signal 
of an epoxy piece used to position the sample). In the 
detail of Figs. i: c) and (d), corresponding to an applied 
field of 4 T, it may be appreciated a rounding extending 
from Tc ~ 28.1 K to ~33 K. This interval is well be¬ 
yond the transition half-width (^0.4 K), so the rounding 
may be attributed to superconducting fluctuations. This 
effect prevents determining the superconducting parame¬ 
ters as in the case of the resistivity, and fluctuation effects 
will be taken into account to analyze the magnetization 
data. 

The fluctuation contribution to the magnetic moment 
was determined through 

mfi{T) = m{T) - mB{T), (5) 

where msiT) is the background contribution coming 
from the samples normal state and to some extent from 
the sample holder. It was obtained by fitting a Curie-like 
function 


mB{T) = A + BT+^ ( 6 ) 

to the raw data in a region between 35 K up to 42 K 
{A, B and C are free parameters). The upper limit was 
chosen to avoid a small upturn (in the scale of 10“^ emu) 
observed above ^45 K in all the measurements, and that 
may be attributed to an unavoidable amount of oxygen 
in the sample space. This upturn cannot be removed by 
successively pumping and venting with helium the sample 


space. The resulting fitting parameters followed a linear 
Fl-dependence within their uncertainty. 

The resulting mB{T) contributions are presented as 
solid lines in Fig. [51 The dashed lines represent the uncer¬ 
tainty associated to the fitting procedure. The resulting 
fluctuation magnetic susceptibility, Xfl(T) = mfi(T)/HV 
(where V is the sample volume), is presented in Fig. [T] 
The uncertainties associated to the background are only 
~ 1% (for H II c) and ^ 5% (for H || ab) in the lower 
bound of the Gaussian region (solid data points in Fig.[7|), 
and negligible in the critical region. Although, as ex¬ 
pected, they are significant for temperatures near Tonset: 
~ 30% for H II c, and ^ 100% for H || ab at 30 K. Notice¬ 
ably, Xfl is anisotropic, being significantly larger in am¬ 
plitude when H A ab. Also, the xn amplitude decreases 
with H, which is an indicative that the fields used in 
the experiments are large enough as to enter in the finite 
field (or Prange) fluctuation regimei^ The quantitative 
analysis of the data would then require using theoreti¬ 
cal approaches valid beyond the zero-field (or Schmidt) 
limit. 


V. ANALYSIS OF FLUCTUATION EFFECTS IN 
THE MAGNETIZATION 


A. Gaussian region above Tc{H) 


It has been shown that the presence of several bands 
contributing to the superconductivity in these materials 
may affect the fluctuation-induced observables through a 
renormalization of the coherence length amplitudes, but 
without appreciably affecting their functional form with 
respect to the single band case<^ Thus, our Xfl(T, iJ) 
data will be analyzed in terms of a GL approach for 
single-band three-dimensional anisotropic superconduc¬ 
tors (3D-aGL approach). The possible presence of multi¬ 
band effects will be probed through differences between 
the resulting coherence lengths and the ones previously 
determined from the field dependence of the resistive 
transition. 

In terms of the 3D-aGL approach the fluctuation mag¬ 
netization Mfl of an anisotropic superconductor in pres¬ 
ence of a finite applied magnetic field is given by^“— 


M^iT,H) = - 


-InF 


kBTj 

T^4>O^abi0) Jo 


dq 


c — e 
2h 


InF 



(7) 


for H T ab, and 


Ml{T,H) = -Mjr{T,H/x). (8) 

7 


for H II ab. Here F and tp are, respectively, the gamma 
and digamma functions, e = ln(r /Tc) the reduced tern- 
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FIG. 6. a,b) Overview of the temperature dependence of the magnetic moment above Tc for all magnetic field amplitudes 
and orientations investigated. The solid lines are the background contributions obtained by fitting a Curie-like function above 
35 K. c,d) Detail around Tc corresponding to the measurement with = 4 T. The shaded areas represent the fluctuation 
effects, extending up to 5 K above Tc and with an amplitude clearly anisotropic. The dashed lines indicate the uncertainty in 
the background contribution. 


perature, h = iJ/[(^o/27r/Xo^^f,(0)] the reduced magnetic 
field, and c is a total-energy cutoff constant The cutoff 
was introduced to take into account short-wavelength ef¬ 
fects, which may be relevant in particular at high reduced 
magnetic fields or temperaturesi^ In view of Eq. o, 
c equals the reduced temperature at which fluctuation 
effects appear, Tonset ~ 33 K. So, in subsequent anal¬ 
yses we will use c = ln(Tonset/7c) ss 0.16, a value 
slightly smaller than the one found in other Fe-based 
superconductors 1 ^5,48,51 

Eqs. 0 and (l 8 |) are expected to be applicable down to 
the critical fluctuation region. This last is bounded by 
the so-called i/-dependent Ginzburg criterion, which for 
3D anisotropic superconductors may be written aS)^“— 




AirkB^.oH 

Ac^c(0)(/)o 


1 2/3 


(9) 


for H T ab, and 


T'^iH)^THH)±T, 


_Ac^c(0)7(/<o. 


2/3 


( 10 ) 


for H II ab, where Ac is the specific-heat jump at T^. 


Let us finally note that Eqs. © and 0 were already 
successfully used to explain the susceptibility round¬ 
ing above of iron pnictides like Bai_xKa;Fe 2 As 2 and 
Ba(Fei_xNia;) 2 As 2 ,— and of other 3D anisotropic su¬ 
perconductors like MgB 2 and NbSe 2 J^‘^ In turn, the 2D 
and 2D—3D (Lawrence-Doniach) versions accounted for 
the behavior of different high-Tc cuprates 

Eqs. (I?!) and ([5|) were fitted to the set of Mfi/H data for 
both H J- ab and H || ab with only two free parameters, 
Cat)(0) and 7 . The best fit is represented as solid lines in 
Fig. [T] The htting interval range from Tonset to the solid 
data points. If the htting interval is extended to lower 
temperatures, the flt quality is considerably worsened. 
In fact, these solid data points are already close to the 
upper bound of the critical fluctuation region. This is 
better seen in the H — T phase diagrams in the insets 
of Fig. [3 where the limit of applicability of the GGL 
approach is indicated by circles, and the shaded areas 
represent the critical regions evaluated with Eqs. m and 
(ITUl) by using AcjTc 0.1 J/molK^^ and the ^c(O) and 
7 values in Table HIl 

The values for the fitting parameters, ^a6(0) and 7 , are 
compiled in Table El together with the ^c(O) value ob- 
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FIG. 7. Temperature dependence of fluctuation-induced mag¬ 
netic susceptibility in the Gaussian region for both H L ab 
(a) and H || ab (b). The lines are the best fit of Eqs. ([Tl and 
respectively, by using 5a6(0) and 7 as free parameters. 
The solid data points are the lower bound of the fitting re¬ 
gion. Insets; H — T phase diagrams for both field orientations. 
The shaded areas represent the critical regions, as obtained 
from Eqs. m and m- The circles (squares) represent the 
lower bound for the applicability of the Gaussian (critical) 
GL approach. 


tained as ^af)(0)/7. The indicated uncertainties account 
for the range in which the fitting parameters may be 
changed without appreciably worsen the fitting quality. 
The ^ab(O) value is in excellent agreement with the one 
determined in Section III from resistivity measurements, 
but the ^c(O) value is a 30 % larger. Such a difference 
is beyond the experimental uncertainty, and will be dis¬ 
cussed in Sec. V.C. 


B. Fluctuation diamagnetism in the critical region 
aronnd Tc{H) 

As a check of consistency of the above analysis we stud¬ 
ied the data in the critical region around the Hc 2 (T) line, 
where the Gaussian approximation breaks downi^ In this 
region, the 3D-GL approach in the lowest-Landau-level 
approximation predicts that and m||(T, iJ) 

follow a scaling behavior, = /_L,||(t_L,||), the scaling 



t (10 ® A'"'® m"'® K^'®) 


FIG. 8. 3D-GL scaling of the magnetization in the critical 
region for H ab (a) and H \\ ab (b). The line in (a) is the 
experimental scaling function for H T ab, while the one in (b) 
for H II ab was calculated from the former through Eq. (nsj. 


variables being^IiZ^ 


Ml 


rn±,\\ = 




( 11 ) 


and 


1-L,II = 


T-Tc ’"(i?) 


( 12 ) 


As always, the indexes T and || stand for H J- ab and 
H II ab, respectively. The scaling functions for H J. ab 
and H \\ ab are related through^ 


/||( 1 ||) = 




v5/3 


(13) 


By assuming a linear temperature dependence of the 
upper critical fields, the scalings for both H J. ab and 
H II ab depend on only three parameters: T^, and the 
H^(T) and slopes. The best scalings, presented 

in Fig. [HI were obtained with the values for these param¬ 
eters presented in the last row of Table HIl The indicated 
uncertainties represent the range of values for which the 
scalings are not appreciable worsened. The correspond¬ 
ing Caf)(0) and 5c(0) values, obtained by using Eqs. ([ 21 - 
El), are also compiled in Table m In spite of the larger 
uncertainties, these values are consistent with the ones 
found in the analysis of the Gaussian region. Other in¬ 
dications of the consistence of the present analysis are: 
i) The scaling is acceptable down to t||__L ^ —3 x 10“®, 
a value in good agreement with the lower bound of the 
critical region (it leads to the squares in the H — T phase 
diagrams in the insets of Fig. Ej). ii) The experimental 
scaling functions (solid lines in Fig. [5]) are related to each 
other as predicted by Eq. m, see the caption in Fig. [5] 



























TABLE III. Summary of the superconducting parameters in 
the literature for BaFe 2 (Asi_a;Pa ;)2 with x near the optimal 
(maximum-Tc) value. 


Tc ModHci/dT ^JLodHljdT Cab(O) 
(K) (T/K) (T/K) (nm) 

5c (0) 
(nm) 

observable 

Ref. 

30.6 

-1.72 

-3.66 

2.50 

1.17 

p(T)h 

79 

30.0 

-2.1 

- 

2.30 

- 

C{T)h 


29.0 

-2.1 

- 

2.30 

- 

C{T)h 


28.1 

-2.23 

-5.73 

2.29 

0.9 

C{T)h 


28.4 

-2.1 

- 

2.33 

- 

C{T)h 


30.5 

-1.67 

-2.41 

2.54 

1.75 

LCR 



iii) The smooth temperature dependence of the magne¬ 
tization around Tc{H) (i.e., around t = 0), see Fig. 8 , 
justifies the seemingly negligible role of the possible T^- 
inhomogeneities in the sample. 

C. Discussion of the results 

As it is shown in Table im the ^a 6 ( 0 ) values deter¬ 
mined from the resistivity measurements are in excel¬ 
lent agreement with the one derived from the fluctua¬ 
tion magnetization. However, there is a significant differ¬ 
ence between the corresponding ^c(O) values. This effect 
cannot be attributed to an incomplete superconducting 
volume fraction in the stack used for the magnetization 
measurements, because ^ab{0) and 7 (or ^c(O)) affect the 
amplitudes of and Mjj, but also their dependence 
on H. The ^o&(0) and ^c(O) values in the literature for 
BaFe 2 (Asi_a;P 3,)2 with values near the maximum one, 
are summarized in Table m These data were obtained 
from the Tc{H) dependence, as probed by dc and ac elec¬ 
trical transport and by a static property like the specific 
heat. The ^at)(0) values in these works are close to each 
other and to the values obtained here (see Table |ll|. The 
^c(O) value obtained from dc resistivitj^ is also in excel¬ 
lent agreement with our data obtained from the same ob¬ 
servable. However, a notable difference is found with the 
^c(O) values derived from other observables (0.9-1.75 nm). 
These results may be somewhat affected by the uncer¬ 
tainty associated to the use of a criterion to determine 
Tc{H), but could indicate a possible dependence of Cc( 0 ) 
on the observable used to determine it. 

Our present results may be due to the presence of sev¬ 
eral bands contributing to the superconductivity in the 
material under studyIn fact, a recently proposed 
two-band model predicts that the functional forms of 
dominating divergences of the fluctuation-induced spe¬ 
cific heat and electrical conductivity do not change with 
respect to the single-band case, but their relative am¬ 
plitudes are affected mainly by a renormalization of the 
c-axis coherence lengthi^ It has been suggested that this 
effect could be behind a dependence of the upper critical 
fields of FeSeo.sTeo.s on the observable used to determine 



doH (T) 


FIG. 9. Anomalous H-dependence of the fluctuation magne¬ 
tization for temperatures just above Tc and for both H ah 
(a) and H \\ ab (b). As a reference, the dashed lines represent 
the ideal diamagnetic response (evaluated taking into account 
the demagnetizing factors of the sample under study). The 
anomalous upturn appears at temperatures within the transi¬ 
tion width (0.8 K full width). Inset: Temperature dependence 
of the upturn magnetic field, Hup. 


them. In particular, the Hc 2 values derived from mea¬ 
surements of the specific heat are systematically larger 
than the ones obtained from the electrical resistivity, and 
these last are in turn larger than the ones obtained from 
magnetic torquei^^— It would be then very interesting 
to check whether an expression for Mfl(T, iJ) in a two- 
band superconductor (still unavailable) would also quan¬ 
titatively account for the differences found here with the 
data obtained from resistivity measurements. 


D. Enhanced diamagnetism in the low-fleld region 
just above Tc 

For completeness, we have also explored this region of 
the phase diagram. In contrast with the linear behavior 
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predicted by Gaussian GL approaches under low fields 
{h e), the Mf[{H) isotherms just above present 
an anomalous upturn at a field ^oH^p ~ 10“^ T for 
both field orientations (see Fig. IH]). Below H^p the mag¬ 
netic susceptibility is orders of magnitude larger than the 
one associated to Gaussian superconducting fluctuations. 
Above Hup the Mg amplitude decreases and, consistently 
with Fig. [71 for fields about 1 T the conventional GGL 
behavior is recovered. 


A qualitatively similar anomalous behavior has been 
observed in high-Tc cuprates like La 2 -xSrxCu 04 ^“— 
and Yi_a;Gaa;Ba 2 Gu 30 y and more recently in Fe-based 
superconductors like SmFeAsOi-xFa, (Smllll)^ and 
Ba(Fei_a;Rh 2 ;) 2 As 2 (Rhl22)^. The similarities between 
our results and these previous works in iron pnictides are 
even quantitative. For instance, the isotherms measured 
at temperatures near Tc present very similar H^p and 
Mup values (in the case of Smllll one has to take into 
account its 2D nature and that the sample is granular). 
So, to a large extent the effect in BaFe 2 (Asi_a;Px )2 has 
the same origin than in Rhl22 and Smllll. In Refs. 1^ 
and m it has been proposed that, like in the HTSG, the 
anomaly is due to the presence of important phase fluc¬ 
tuations. Alternatively, it has been proposed that it is a 
consequence of a Tc distribution, which in cuprate or Fe- 
based superconductors would be associated to their non- 
stoichiometric nature.— In the case of BaFe 2 (Asi_a;Pa ;)2 
some indications support this last explanation; On the 
one side, below H^p the magnetic susceptibility is a sig¬ 
nificant fraction of the perfect diamagnetism (indicated 
as a dashed line in Fig.lH]), as if a significant volume frac¬ 
tion of the sample were actually in the Meissner region. 
On the other, the effect is limited to a range of tem¬ 
peratures above Tc of the order of the transition width. 
Finally, the H^p value at the average Tc is consistent with 
the lower critical field of the highest-Tc domains, which 
may be approximated by —{dHci/dT)T^ATc, and that is 
in the 10“^ T range— 


An argument supporting the possible presence of phase 
fluctuations in Smllll and Rhl22 is that H^p increases 
with the temperature above Tc, while the lower critical 
field decreases . Such an Hup increase with T is 
not observed in BaFe 2 (Asi_xPx )2 (see inset in Fig. |5|), 
in spite that the transition width of our sample is com¬ 
parable or even smaller (in relative terms) than those in 
Refs, m and [13. As a consequence, we may assume that 
phase fluctuations do not play a role in this compound, 
or at least they are much less relevant than in Smllll 
and Rhl22. 


VI. CONCLUSIONS 

We have presented measurements of the magnetiza¬ 
tion and electrical resistivity around the superconduct¬ 
ing transition of the isovalently-substituted iron pnic- 
tide BaFe 2 (Asi_a;Pa :)2 with x « 0.35, near the optimal 
(maximum-Tc) value. The measurements were performed 
with magnetic fields up to 9 T in the conductivity and 
6 T in the magnetization, in both cases applied in the 
two main crystal directions (parallel and perpendicular 
to the ab layers). Due to the relatively large normal- 
state electrical conductivity near Tc, the fluctuation ef¬ 
fects in this observable are almost unobservable. As a 
consequence, the resistive transitions are very well de¬ 
fined, allowing to determine with accuracy the tempera¬ 
ture dependence of the parallel and perpendicular upper 
critical fields, and the in-plane ^ab(O) and out-of-plane 
^c(O) coherence lengths. On the contrary, the magneti¬ 
zation presents an appreciable diamagnetic contribution 
just above Tc that may be attributed to superconduct¬ 
ing fluctuations. These data were consistently analyzed 
in the different fluctuation regimes (Gaussian and criti¬ 
cal) in terms of existing Ginzburg-Landau approaches for 
single-band superconductors. While the resulting ^ab(O) 
is in excellent agreement with the one obtained from the 
field dependence of the resistive transition, ^c(O) is sig¬ 
nificantly larger (^ 30%). A dependence of the coher¬ 
ence lengths with the observable used to determine them 
was also found in other Fe-based superconductors— “— 
Recent theoretical approaches show that the presence of 
several bands contributing to the superconductivity af¬ 
fects the fluctuation effects mainly through a renormal¬ 
ization of the c-axis coherence length— Thus, we sug¬ 
gest that our present results could be a consequence of 
the multiband nature of the compound under study. An 
expression for the fluctuation magnetization in two-band 
superconductors would be desirable to confirm this pro¬ 
posal. Finally, from the analysis of the anomalous up¬ 
turn in the M{H) isotherms just above Tg, we find that 
phase fluctuations may play a negligible role in this com¬ 
pound, contrary to what is proposed for other Fe-based 

superconductors— 
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